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We propose a kinetic theory of transport in metals embedded with ferromagnetic nanoclusters, 
describing the dependence of the form of magnetoresistance on anisotropy characteristics of the 
ensemble of clusters and the influence of the electron spin depolarisation by clusters. We note that 
this effect is strongest when all the clusters have the same intrinsic easy-axis anisotropy. 
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. Numerous applications of magnetoresistance (MR) phenomenon fuel a continuous search for new materials, including 
' various hybrid systems [ij. In this pubHcation, we propose a kinetic theory of high temperature MR in a normal metal 
embedded with ferromagnetic metallic nano-clusters (FmnC) 0- 
H [ The presented below analysis extends the earlier phenomenological study of MR effect due to FmnC 01 developed 
. in application to granular normal-ferromagnetic metal films 01 • Following Zhang we study clusters using classical 
^5 I spin approximation. We also assume that both the magnetisation and the electron exchange spin-splitting J are 
, homogeneous within each cluster (due to their small size, typically 10^ — lO^^ atoms per cluster 0, IH 0) and we 
psj characterise the ensemble of FmnC using the angular distribution of their magnetic moments m = fil, which depends 
on intrinsic anisotropy, external magnetic field and temperature. This study is restricted to materials with such 
separation between clusters that one can neglect correlations between them and assume that all clusters are of the 
same radius, rd. 

Due to band mismatch between underlying materials, each FmnC would generate a scalar potential well and have 
a different electron mass from that in the normal metal, in addition to a spin dependent (exchange) perturbation. An 
interplay between all above-mentioned factors makes electronic scattering from each individual FmnC spin-dependent 
and the resistance of a composite material sensitive to the mean degree of polarisation of clusters {Iz). In agreement 
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' with earlier publications 0,0,Q|, we find that in the low field regime {fiB/kT <C 1) MR is quadratic in magnetic field, 
^ • A (B) = ~ —B^^ whereas the detailed form of MR is shown below to depend on both mean polarisation of 

FmnCs and its variance. We also analyse the infiuence of magnetic anisotropy on MR, A (B) and find that the effect 
is strongest when all clusters in the ensemble have collinear easy-axis anisotropy. 
To describe the MR in a metal with FmnCs we study the quantum kinetic equation. 



9tp + V • Vp + e (E • v) d,p - ^ [<73, /}] = /o [p] + {I [p]) 



J> . for the electron spin density matrix p — Po+X^i Pi^i approximated by pa (p) — Pa {p) + ^''na (p), where cr = {ai, a2, (T3) 
I are PauH matrices, lub — Beh/mc is the electron spin precession frequency in normal metal and p is the isotropic part 
of the electron distribution in momentum space. All details of the electron spin and phase evolution inside clusters 
are incorporated into the collision integral / [p] , whereas momentum relaxation due to non-magnetic impurities and 
phonons in the normal metal is taken into account by Iq [p] — [p (p) — p (p)] r^^. 
Tij- ' One FmnC produces the following scattering matrix, 

o ■ 

S = a + b{<r-l); a = ^^^,b=^—^, (1) 
• where S^^, and S^^, characterise the angular dependence of scattering matrices of electrons with spin parallel and 



antiparallel to the cluster polarisation 1 [p, p' are electron momenta before and after scattering]. The collision integral 
taking into account a group of collinearly polarised clusters with the concentration n (I) has the form 



Spp'Pp'Sp,p - - {^Spp'Sp,pPp + PpSpp'Sp,^ 
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CZ ' where v is Fermi velocity. This expression incorporates the kinetics of both electron scattering and spin precession, 
in particular since S matrix takes into account spin-dependent forward scattering. Due to the angular distribution of 
cluster polarisations, the ensemble-averaged coUision integral, {I[p]) would also describe spin relaxation of electrons. 
All of the above-mentioned three generic elements of electron spin kinetics would be present even if the distribution of 
cluster magnetic moments is axially symmetric about the direction of an external magnetic field B = Be^- Therefore, 
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in the following analysis we assume that {Ix) = {ly) ~ and the tensor {Ulj) is diagonal with (l"^) — (Z^) — 1 — {if}- 
The resulting ensemble-averaged collision integral reads as 
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where density matrix p was decomposed into a singlet(po) and triplet(pi) parts [here p = p(p), p' = p(p')]i which 
decouples the kinetic equation into two groups. The first group describes the evolution of only diagonal density matrix 
components and involves momentum {t^^) and spin {t~^) relaxation rates, 



dtpo + f V • r?o = 

dtPz + • 773 = - (1 - {ll)) t7o'p3 (3) 

vVpo + evEd^pT (e) = -Tp_^Vo - V- ^'^) 
vVp3 = -T-} {h) rjo - (r-+i + (1 - (ll))) V3, 
V+ = '-^JdOsmehs+f + \S-f] (1 - COS0) + i 

' ' ^ /rf^sin^ ( |5+|' - |^-|') (1 - COS0) , (4) 
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The second group, 
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where C = Pi + and $ = ^71+ 1^72, takes into account an average of electron spin precession in the (xy)-plane with 
frequency 
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Eq.lEJ describes the kinetics of charge (po) and spin (pa) densities in a macroscopic sample. Without any electric 
bias this would be diffusion Q described by equations 

dtpo = -DeV^po - -DaV^ps , , 

dtP3^DsV'p3~DA^'po-{l-{ll))T-o'p3 ' ^''> 
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A homogeneous solution of Eq. IpJ in the presence of electric bias produces the current-field relation for both charge 
and spin components je = crgE and js = cr^E , which determines electric(e) and spin(s) conductivities, 

ae = e^iy (ef) D^, as ^ - {h/2) ev {tp) Da- (7) 
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Table I: Example of estimations of Ac in a spherical cluster model for large clusters prd S> h and strong exchange J 



A magnetic field dependence of average cluster polarisation {Iz) and mean square (Z^) generates a finite magneto- 
resistance (MR) of a composite material, 

A(B) = ^-1. (8) 

The resistance change which can be achieved across an "infinite" magnetic field interval (providing a complete 
uniaxial polarisation of all clusters) is given by 

^max = [Tp+/Tp^f ■ (9) 

The maximal MR effect, /S.max is the largest when both potential and magnetic scattering potentials are of the same 
order, and it is suppressed by additional non-magnetic scattering processes increasing versus t~\ 

The entire MR curve can be described in terms of magnetic field and crystalline anisotropy dependent values of 
(i.) <mtl Ifi). 

When magnetic clusters represent the only source of scattering in the material, the parameter X in En. ltT71|l is 
determined by the cluster size and material composition. For large clusters with radius rd ^ h/pp, where pp is 
the Fermi momentum in the metal all relaxation rates are proportional to the clusters' geometrical cross-section, 
'''a^ ~ nipp'^ (P'^'ci)^ Aa, (ck = p+,p— , sl). Tabled illustrates values of these parameters for various ratios between the 
exchange J and Fermi energy Ep in the ferromagnet and the ratios mF/rriM between band masses in the normal and 
ferromagnetic metals. Here, we use a simplifying model with parabolic dispersion in both materials with equal densities 
of conduction band electrons. The latter constraint determines a spin-independent band mismatch U = Ep — Ep at 
the interface which is taken into account as a scattering potential. The data in Table I indicate that the depolarisation 
parameter X is not large {X < 0.5). For small clusters {prd ^ h), scattering is isotropic, thus in the integral in 
Eq.l^l is suppressed, as compared to the momentum relaxation rate t~^, thus making X ■^1. Addition of potential 
scatterers and a contribution from e-phonon scattering produces a similar effect on the depolarisation parameter. 

Below, we compare the magnetic field dependence A [B) for three ensembles of clusters: (a) intrinsically isotropic 
FmnCs; (b) identical clusters with colinear easy axes ahgned with the external magnetic field; (c) ensemble of clusters 
with randomly oriented easy axes. In all of these three cases we assume the same material composition, size of clusters, 
surrounding medium (i.e. fixed values of rp±, Ts„ and X), and thermal equilibrium of angular distribution of clusters. 
These three situations will be identified with the following form of free energy of magnetic subsystems, 

r -liBU (a) 
F = -n'B-\ + Fn = I -liBU-all {b) (11) 

{ -fiBl, -a {I- lof (c) 

where lo determines the direction of single cluster easy axis. In cases (a) and (b) the values of (Iz) = h and (ll) = 11 
are determined by the thermal average 



A = N-' / dVliAe-T^, N = / df^ie"^, (12) 
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Figure 1: A (B) for different cluster subsystems a) intrinsically isotropic clusters; b) identical clusters with colinear easy axis 
in direction of external magnetic field; c) ensemble of clusters with randomly oriented easy axis; kT = 0.1, a = 1. Solid line 
shows results for X = 0.5 and dashed line for X = 0. 



where T is temperature. For ensemble (c) with isotropic distribution of easy axes, an additional averaging over 
directions of vector Iq in En. (fTTH is needed (1^) — {{h)) and (Z^) = (^(I'^'j'j, where (((. ..)})= / dfli^ (...) /Att. 
For intrinsically isotropic clusters, ensemble type (a), 

{lz)a ^ (xcothx — 1) /x, X = B^/kT, 

{ll)^ = (x^ -2j;cothj; + 2)/j;^ (13) 

The resulting form of MR is illustrated in Fig.l by the curve (a). The influence of intrinsic easy-axis anisotropy of 
clusters on MR depends on the temperature regime and the type of cluster ensemble (b) or (c). At high temperatures 
kT > a, the existence of an anisotropy axis common for all clusters results in a small correction to {Iz), (ll) in 
ensemble (b), 



thus leading to small corrections to A{B). Here L\{x) — [2x^ (coth^ x — l) + 2a::cotha; — 4] /x"^ and L\{x) — 
[4x^ (2 — coth^ x) — IGxcotha; + 20] /x^. In ensemble (c), which has an isotropic distribution of magnetic moments 
at B = 0, such a correction is even smaller and appear only in the second order in a/kT<l. 

At low temperatures, kT < a, A {B) was studied numerically. The obtained results for ensembles (b) and (c) are 
compared in FigQlwith the MR in ensemble (a). FigQ] shows that a pronounced easy-axis anisotropy enhances the 
low-fleld MR effect, whereas an isotropic distribution of easy axes suppresses MR via spin relaxation. The inset to 
Fig.l demonstrates the low-fleld behaviour in all three cases, which complies with the asymptotic expansion of A {B) 
in Eq-Cni in powers of ^iB/kT < 1. Since at small flelds e^^'-/'='^--P'o/fcT ^ ^-Fo/kT _^ ^S;^ 



^ ^ kT Jdfliexpi-Fo/kT) kT^''B=o 



hence the low-fleld MR is quadratic. 



A (B) « — ^ ^ ^niaa; {^1)^=0 



kT) 1 + {1-{11)b^^)X- 

Due to isotropy of FmnC ensembles (a) and (c) at -B = 0, {Idj) b=o ^ ^ b=o^(^'^)c b=o ~ whereas 



The low-field MR in each of those three cases is 



A{B) ffBY \ (9 + 6X) ^ (a,c) 



In summary, we propose the theory of classical magnetoresistance in a normal metal with embedded ferromagnetic 
nanoclusters. Our analysis shows that in a thermodynamic equihbrium the MR effect is largest when all clusters 
have the same intrinsic easy-axis anisotropy. Thus it may be advantageous for applications to engineer anisotropy as 
sketched in FigQ] despite the inconvenience which may be caused by hysteresis in their magnetisation dynamics. 
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